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NOMENCLATURE 

a - Transformed length in the transformed plane for 0.5 H in the phys- 
ical plane 

/ - Fractional volume of burned medium in a cell 

- Side length of a cell 
H - Width of jet exit 
i - VM 

m - Index for any vortex blob 

fi f - Unit vector normal to the flame front 

P - Non-dimensional pressure 

f - Position vector (=x+iy) non-dimensionalized by H 
fj - Non-dimensional position vector of flame front 
t 0 - Blob core radius 

R, - Reynolds number 

s c - Strength of combustion generated volumetric source 

S u - Non-dimensional burning rate 

t - Time non-dimensionalized by H/ U 0 

t x - Time instant for starting combustion 

t 0 - Time instant for starting ignition 

t p - Time instant for peak pressure 

A t - Non-dimensional time increment 

u - X-component velocity in the physical plane 

u - Velocity vector non-dimensionalized by the jet inlet velocity 

U m - X-ccmponent velocity in the transformed plane 

U g - Jet inlet velocity 

v - Y-component velocity in the physical plane 
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v b - Volume of burned medium 
V c - Volume of a cell 

V m - Y-component velocity in the transformed plane 
V u - Volume of unburned medium 
fV b - Complex potential of a blob 

- Complex potential of a blob image 

W c - Complex potential of a combustion generated volumetric source 
w cg - Complex potential of a combustion generated volumetric source 
image 

W m - Complex potential of any vortex blob 
W, - Complex potential of a source 

x - X coordinate in the physical plane 

y - Y coordinate in the physical plane 

Z - Complex position coordinate (=x+iy) in the physical plane 

( - Complex position coordinate in the transformed plane 

<6 - Complex position coordinate of a blob in the transformed plane 

- Complex position coordinate of a combustion generated source 
in the transformed plane 

A - Wave length 

y. - Viscosity 

Pt, - Density of burned medium 

p u - Density of unburned medium 

o - Standard deviation (= Vf» 

u - Vorticity (=Vxtz) 

u b - Vortex strength of a blob 

V - Non-dimensional del operator 

V 8 - Non-dimensional Laplacian operator 
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- Non-dimensional substantial derivative 



CHAPTER 1 


INTRODUCTION 

One of the keys to our energy future is to develop a lean combustion 
engine that can operate stably and economically using lean mixtures and 
producing clean emissions 1 . The concept that a lean combustion engine 
offers distinct advantages , both with respect to pollution as well as 
efficiency, is well known, as attested by the development of stratified 
charge engines. However, efforts to achieve satisfactory performance 
under extra-lean conditions have been so far unsuccessful. Most of the 
difficulties encountered in this connection can be ascribed to an overly 
localized ignition process. That is, the rate of burning the charge depends 
on the flame propagation speed, which in lean mixtures is just too slow 
for satisfactory operation. Thus, attempts to develop a lean combustion 
engine invariably result in problems associated with excessive emissions 
of unburned hydrocarbons due to poor flame propagation. To make 
matters worse, the bulk quenching^ and the wall quenching effects ^ 
enhance the emission index of unburned hydrocarbons. These 
deficiencies can be overcome by the use of multi-point ignition sources 
since high overall rates of burning in a lean charge can be attained by 
having the combustion initiated at a multitude of sites distributed 
throughout the mixtures. The multi-point ignition can be achieved by 
impregnating the charge with active radicals which enhance the induc- 
tion process and speed up the onset of exothermic reactions. This can be 
accomplished by the use of combustion jet ignition system which delivers 
and distributes the free radicals throughout the combustion chamber 
with turbulent jets. 
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The combustion jet ignition system is a promising means for 
enhanced ignition in I. C. engines with lean premixed charge^. The sys- 
tem generates jets of combustion products containing free radicals and 
discharges them into the combustible medium as ignition sources. These 

g 

initial concentrations of free radicals serve as substitutes for high initial 
temperature in the ignition of combustible mixtures. The use of jets of 
free radicals as ignition sources for combustible mixtures offers the fol- 
lowing advantages : 

(a) a controllable depth of penetration, so that combustion can be 
started at any desired location within the charge, 

(b) zonal pre-turbulization enhancing the combustion process, and 

(c) a wide dispersion resulting in multi-point ignition. 

Of particular interest to us is the fact that these three properties of 
combustion jet ignition can improve the thermal efficiency of the 
combustion processes in I. C. engines by promoting high overall rates of 
burning. To illustrate this improvement, consider the following example. 
In an idealized Otto engine, as shown in Figure 1, the constant volume 
combustion process 2-3 is where ignition and inflammation processes 
take place. Prior to the ignition process are the intake process 6-1 and 
the piston compression process 1-2. The intake valve can induce a 
significant amount of turbulence which persists until ignition. The pis- 
ton motion during the compression stroke also can induce turbulence, 
manifested primarily by the roll-up vortices^ These induced tur- 
bulences can strongly influence the ignition processes since they exist 
prior to ignition. After the inflammation process are the piston expan- 
sion process 3-4 and the exhaust processes 4-5 and 5-6. The higher the 
pressure of state 3, the higher will be the thermal efficiency of the Otto 
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cycle. The higher pressure of state 3 can be achieved by using the 
combustion jet ignition which promotes high overall rates of burning. 
Figure 2 shows the comparison of the pressure-time traces obtained in 
closed bomb experiments for conventional spark ignition, and combus- 
t i jet ignition. The third curve in Figure 2 shows theoretically calcu- 
lated pressure-time trace for constant volume combustion using CEC-76 
computer program . The induction period (t^-t Q ) is shorter for combus- 
tion jet ignition than for conventional spark ignition. In addition, the 

combustion duration (t -tj is also shorter for combustion jet ignition 

P * 

than for conventional spark ignition. With shorter induction period and 
combustion duration, combustion jet ignition has higher peak pressure 
due to fast burning and less energy loss. 

A combustion jet ignition system consists of the following five distinct 
processes: 

1. The generation of combustion products containing free radicals. 

2. The delivery and distribution of the free radicals by turbulent 
jets. 

3. The induction period. 

4. The ignition processes. 

5. The inflammation processes. 

The understanding of these five processes are crucial to the develop- 
ments of the combustion jet ignition system, which is one of the proper 
means of controlled combustion to govern the action of free radicals. 
These five processes cover the studies of the mechanism of the genera- 
tion, the delivery and distribution of free radicals, and of the role they 
play in controlling the combustion processes. Reported here are parts of 
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the studies: 

1. Experimental study of the formation processes of the combus- 
tion jets of free radicals. A specially designed igniter with 
prechamber was used to build up and control the stagnation 
pressure upstream of the orifice. 

2. Theoretical study of the formation processes of the turbulent 
jets by using Chorin’s Random Vortex Method (RVM). 

3. Theoretical study of the inflammation processes in the turbulent 
jets by using RVM implemented with the flame propagation algo- 
rithm. 



CHAPTER 2 


LITERATURE SURVEY 

For the study of combustion jet ignition, surveys on bot:* ignition stu- 
dies and turbulent jets were performed. In this chapter, the literature 
survey on ignition studies will be presented first, followed by the litera- 
ture survey on turbulent jets. 

2. 1 Ignition Studies 

Ignition is the initiation of a self-sustained exothermic process of 

13 

combustion. It is different from inflammation which is the initiation 
and propagation of a self-sustained flame. Ignition is an essential 
ingredient of controlled combustion. In order to understand the mechan- 
ism of the ignition processes, several ignition systems are being actively 
investigated by researchers. Among those ignition mt.<.nods, the following 
six were selected for literature survey here: 

1 Shock Induced Ignition Method. 

2. Conventional Spark Ignition Method. 

3. Photochemical Ignition Method 

4. Flame Jet Ignition Method. 

5. Plasma Jet Ignition Method. 

6. Combustion Jet Ignition Method. 

2.1.1 Shock Induced Ignition Method 

The reflected shock tube technique is widely used to obtain data for 
chemical reaction rates in the high temperature range. This te ‘hnique 
provides a test environment of uniform temperature and pressure behind 
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reflected shock waves for studying the ignition processes. Many interest- 
ing and fundamental physical concepts of ignition processes have been 
reported from shock induced ignition studies 14 ' 17 . These physical con- 
cepts of ignition processes are: 

1. Ignition has two modes : (i) a constant pressure process 
representing the case of mild ignitio- and (ii) an adiabatic pro- 
cess representing the case of strong ignition. 

2. The demarkation line between these two modes on the pressure- 
temperature plane of initial conditions is called the strong igni- 
tion limit and is best correlated with the line ( - ■ rJ ) p - constant 

OT 

for hydrocarbon-oxygen mixtures, where r is induction period, T 
is temperature and P is pressure. 

3. Higher energy release rates can generate blast waves within a 
smaller volume of gas in the course of ignition. 

4. Higher temperature sensitivity of the induction time can cause a 
smaller fraction of the gas to be auto-ignited. 

5. The dynamic properties of exothermic processes in combustion 
are due to the evolution of mechanical energy that is manifested 
by pressure waves which affect the flow field. 

2. 1.2 Conventional Spark Ignition Method 

The conventional spark ignition system supplies energy between 10- 
50 mJ per electric spark event of a duration of approximately 1 ms and 
ignites homogeneous and stoichiometric or slightly rich air-fuel mixtures 
without difficulty. The processes of spark ignition depend on many 
parameters such as arc energy, peak voltage, duration of discharge, 
geometry of the spark gap, and the location of the spark relative to the 
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particular geometry of the compressed charge . generally, lean mixtures 

require additional precautions, such as increased electrode distances in 

the spark plug, prolonged duration, higher energy, etc. In a current 

1 fl-20 

detailed study of the new aspects on spark ignition , some interest- 
ing findings of spark ignition were summarized as follows: 

1 . The detailed structure of the spark event was divided into spark 
head and O low phase - , and separated by a transition region. 

2. About 10% of the total energy is used to produce the spark 
plasma within 10 fjs during the spark head phase for every spark 
event. 

3. The power input in the spark head phase is very effective and of 
high importance in the ignition processes. 

4. The temperature of the spark head phase is about 60000 K. This 
high temperature causes the initiation of combustion associated 
with the generation of a pressure wave in the thin spark channel 
(» 3 mm 3 ) which starts to expand with supersonic velocity. 

5. The temperature of the glow phase of the spark event seems to 
occur in a temperature range somewhat higher than 1000 K. The 
energy losses to the electrodes during this relatively long time (» 
1 ms) become more important. Therefore the glow phase of the 
spark cannot be as efficient as the spark head phase with regard 
to energy transferred to the gases. 

6. The starting diameter of the spark channel for ignition is d 1 er- 
mined by the breakdown mechanism. It is possible to increase 
this starting diameter (or starting volume) by temporal redistri- 
bution of the energy of the Capacitor-Discharge ignition without 
changing the total energy supplied to the spark. 
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2. 1.3 Photochemical Ignition Method 

In the studies of chemical kinetics, photochemical techniques were 

developed to experimentally investigate the photolysis processes. By 

21-23 

applying this technique to the studies of ignition processes , some 
observations were obtained as listed in the following: 

1. Photochemical ignition does not impose a temperature rise but, 
instead, increases the concentrations of important reaction 
intermediates and shifts the thermodynamic state of the mix- 
tures into the rapid reaction region via photodissociation. 

2. The ignition and flame stabilization can be achieved by using 

-5 -5 

vacuum ultraviolet radiation (1.45 xlO <A< 1.80x10 cm) 

3. The critical concentration of oxygen atoms for igniting a 
hydrocarbon-air mixture is of an order of 10^ atoms/cm^. 
Higher initial oxygen atom concentrations can reduce the induc- 
tion period and enhance flame propagation. 

4. At input energies of interest to internal combustion engine sys- 
tems and for plasmas generated in certain observed fuels, down- 
stream effects of vacuum ultraviolet radiation make no 
significant contribution to large improvements in the flame ini- 
tiation observed during ignition by hydrocarbon plasmas. 

In addition to vacuum ultraviolet radiation, laser light which is in the 
visible and infrared regions of the spectrum was also successfully used to 
ignite the hydrocarbon-air mixtures. The laser beam is focused at a 
small point within the mixtures to achieve breakdown, thereby generat- 
ing a plasma kernel, which is approximately 2 cm long and 0.5 cm in 
diameter, that acts as a site for the source of ignition. The investiga- 
tions, carried out using an ASTM-CFR engine^ demonstrated the 
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following results: 

1. The laser ignition system has a more rapid pressure rise than 
the conventional spark ignition system for all air-to-fuel ratios 
tested. 

2. The conventional spark ignition system would not run with an 
air-to-fuel ratio above 22.5:1, while with the laser system the 
maximum air-to-fuel ratio obtained was 27.8:1. 

3. Laser systems can reduce specific fuel consumption and have 
less HC and CO emissions, but have more N0 X emissions. 

2. 1 .4 Flame Jet Ignition Method 

The flame jet ignition system can generate a jet of flames that 
penetrate deeply into the main charge to initiate the inflammation 
processes. The system initiates the combustion in the prechamber first 
and causes a pressure difference between the prechamber and main 
chamber which are connected by a channel. Before the completion of the 
combustion in the prechamber, the flames are injectf d, due to the pres- 
sure difference, into the main chamber through the inlet channel con- 
nected to the prechamber. Depending on the pressure difference, the 
diameter and the length of the inlet channel wall, the jet of flames can be 
either laminar or turbulent. As the jet of flames exits the inlet channel at 
a relatively high velocity, bulk quenching can occur. Furthermore, the 
jet of flames experiences shearing due to the effects of the channel walls. 
As a consequence, the jet of flames can be a jet containing a large 
number of small size turbulent flamelets or bulk quenched turbulent ker- 
nels if a high enough degree of turbulence can be achieved. 
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The flame jet ignition system distributes the jet of flames, which can 
be bulk quenched, throughout the combustion chamber to initiate the 
inflammation processes. The inflammation processes initiated by this 

C neon 

system have been extensively studied by researchers 0, co ", but the 
results of these studies are still controversial and need further investiga- 
tion. Gussak et a L define this inflammation process initiated by 

flame jets as the "LAG" process, i.e. avalanche activation of combustion 
process. Applying his so-called "LAG" process to prechamber engines, 
Gussak claimed the following: 

1. The "LAG” process has nothing to do with the frontal mechanism 
of combustion. It is associated with the mechanism of bulk 
inflammation of a combustible mixture. It has an inflammtion 
delay period. 

2. Though the temperature of flame jets is lower than that in the 
electric discharge of a spark ignition, nevertheless, flame jets 
can reliably ignite the mixture in the main chamber. 

3. Flame jets can ignite a wider range of lean mixtures in the main 
chamber by a factor of 1. 5-2.0 than a conventional electric spark 
can. 

4. The main chemically active center in the combustion products of 
a rich methane-air mixture is atomic hydrogen. 

5. A steady flame jet can stabilize and promote combustion of the 
main flow of a combustible medium in a combustion channel with 
no stabilizer. Ti;is steady flow of flame jets impinges perpendicu- 
lar to the main flow at a small subcritical pressure difference 
(about half the critical value) between the prechamber and the 
main chamber. 


v„ 
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6 For an optimal design of a prechamber system, the prechamber 

has a relative volume of 2-3% of the main chamber volume and is 

equipped with outlet channels whose relative cross section is 
2 3 

0.03-0.04 cm per cm prechamber volume. Under these 
optimal conditions, a mixture of equivalence ratio 1. 4-2.5 in the 
prechamber can provide the highest specific power of the 
prechamber engine, minimum specific fuel consumption, max- 
imum stability and the widest range of effective leaning of mix- 
ture in the main chamber. 

Besides Gussak et a L, Wang et ai in China have also studied flame 
jet ignition with an application to a simulated rotary combustion engine. 
They made the following statements: 

1 . By using lean mixtures, the flame jet igniter produces jets of 
flame which either ignite the mixture stream in the main 
chamber or extinguish at the mouth of the prechamber. 

2. By us 4 rich mixtures, the flame jet igniter produces jets of 
flame which become flamelets, traveling along the jet path fine 
for some distance from the outlet. They either extinguish or flash 
to ignite the mixture stream in the main chamber suddenly. 

3. Under t\ : same stream velocity in the main chamber, flame jet 
ig T ”t.jn using lean mixture is due to high temperature effect, but 
flame jet ignition using rich mixture is due to active radicals 
effect. 

Based on their theoretical chemical kinetic calculations of the igni- 

e pq no 

tion mechanism, Creighton et at < ’ do not agree Gussak’s theory for 


the "I ;.G" process engine. They claimed: 
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1 The flame propagation is dominated by diffusion of radicals which 
eliminates any induction period. 

2. The buildup of radical concentrations governs the induction 
period, and the rate of fuel consumption is determined by the 
steady state radical concentrations. 

3. The ignition temperature can be decreased by addition of radi- 
cals. This effect becomes especially noticeable when the initial 
CHg concentration is at a level exceeding 1% of the initial con- 
centration of fuel. 

4 The calculations show lower radical concentrations for the "LAG" 
process than those Gussak claimed, but the results are of the 
same order of magnitude. 

5. The "LAG" process is a flame propagation process which does not 
depend on induction type chemistry. The flame produces enough 
radicals to achieve steady state radical concentrations ahead of 
flame. 

Finally, Furukawa et at?' 1 in Japan investigated the controllability of 
the time required for combustion in the main chamber using turbulent 
flame jets. The results of their studies show that the flame jets first enter 
the main chamber and then start inflammation. As the orifice diameter 
decreases, the penetration speed of the leading front of the flame jets 
becomes higher. Moreover, the combustion duration in the main 
chamber tends to decrease rapidly for the rich mixture of propane-air 
and n-butane-air than others. 
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2. 1.5 Plasma Jet Ignition Method 

The plasma jet ignition method produces a high temperature plasma 
in a confined, recessed cavity provided with an orifice. The high tempera- 
ture plasma is generated so rapidly that the cavity is pressurized, caus- 
ing a jet of plasma to be discharged through the orifice into the combus- 
tible medium as an ignition source. The electric power supply for plasma 
jet ignition system discharges at a relatively low voltage and high current 
through the spark which is generated in a conventional manner by a high 
voltage, low current ignition system. 

Considerable research work^*^ has been done on the studies of 
plasma jets and their ignition characteristics. Some observations were 
obtained as listed in the following: 

1. Plasma jets enter the main chamber in the form of a turbulent 
plume which is imbedded in a blast wave headed by a hemispher- 
ical shock wave. 

2. The blast wave effects are dissipated by the time combustion 
starts. After a delay typically of an order of 1 msec, ignition 
takes place in the turbulent zone of the plume. 

3. The depth of penetration of the jets is a function of their initial 
velocity. 

4. The burning speed is initially quite high and decreases monotoni- 
cally as the flame kernel expands. 

5. With provisions made to fill the cavity with different feedstocks, 
the most effective for ignition were hydrocarbons, among which 
those initially in liquid state were particularly effective. 
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Plasma jets can ignite gaseous mixtures below the normal 
flammability limit. 

7. A plasma jet ignition system requires relatively high electric 
power, eg. 1 J. 

8. The HC and NO emissions are increased and the CO emission is 
slightly decreased. 

2. 1.6 Combustion Jel Ignition Method 

The combustion jet ignition system can generate jets of combustion 
products containing free radicals and discharge them into the combusti- 
ble medium as ignition sources. 

The combustion jet ignition method is very similar to the flame jet 
ignition method. The only difference between these two methods is that 
the latter method generates a jet of flames and the former method pro- 
duces a jet of combustion products. The compositions of flame jets are 
incomplete combustion products containing flames, but the compostions 
of combustion jets are complete combustion products. In order to study 
the flame jet ignition mechanism, one has to understand the effect of the 
degree of incompleteness of combustion products on ignition processes. 
In doing so, a study of combustion jet ignition, which uses a jet of com- 
pletely combusted products, has to be performed first. Then the degree 
of incompleteness for jets of combustion products is varied to study 
flame jet ignition. 

The combustion jet ignition method is also very similar to the EGR 
ignition method which distributes a virtually infinite number of ignition 
sources throughout the fresh charge to initiate combustion in bulk at a 
proper rate of burning. The EGR ignition system is a system using a 
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sufficiently large amount of exhaust gases to ignite their mixture with the 
fresh charge. Such a system has been developed, up to commerical 
application, by Onishi at aZ. 4 ^ with a remarkable proof of clean emissions 
and higher thermal efficiency. 

Up to now there is no literature showing the characteristics of the 
combustion jet ignition method. But, with the knowledge obtained from 
both flame jet ignition and EGR ignition studies, a special combustion jet 
igniter with a prechamber for generating combustion jets has been 
developed at University of California at Berkeley. The igniter can also 
generate jets of flames with varying degrees of incomplete combustion 
for studying flame jet ignition. 

2.2 Turbulent Jets 

Turbulent jets are the major means- of controlled combustion for the 
delivery and distribution of free radicals, which are essential intermedi- 
ate constituents. They serve as multi-point ignition sources to improve 
the thermal efficiency of the combustion processes in a lean combustible 
mixture. In order to understand the fluid mechanic properties of tur- 
bulent jets, a survey of turbulent jets was performed. 

Turbulent jets have been extensively studied 43 " 5 ® both theoretically 
and experimentally. The theoretical studies of plana or axi-symmetric 
steady turbulent jets were performed traditionally using self-similar 
approximations 43 ’ 45 . A recent study 4 ® used the "pressure-correction" 
numerical method to approximate two-dimensional, time-dependent tur- 
bulent jets. For the case of steady turbulent jets with combustion, e g. 
fuel jets 45,4 ®, the plane or axi-symmetric fuel jets were approximated 
using self-similr method implemented with the Schwab-Zeldovich 


transformation technique. All the studies discussed above do not show 
the detailed structure of turbulent jets such as the large-scale eddy 
structure of turbulent jets. 

Turbulent jets have also been investigated experimentally. A variety 
of experimental data on plane and axi-symmetric turbulent jets without 
combustion have demonstrated the self-similar feature of turbulent 
jets . More recent experimental interests are focused on the formation 
of the turbulent plume with and without combustion This experimen- 
tally observed turbulent plume has been actively investigated by many 
researchers as mentioned in sections 2.1.4 and 2.1.5, in an attempt to use 
it as an ignition tool. A detailed experimental study of the axi-symmetric 
turbulent jets with and without combustion has shown the existence of 
the large-scale eddy structure of turbulence^ 7 . The large-scale eddy 
structure of turbulent jets, as shown in Figure 3, is organized in the form 
of a staggered array of vortices of alternating circulation. This large- 
scale eddy structure of turbulent jets, as with other large-scale eddy 
structures in other turbulent shear flows, dominates the mixing and 
entrainment processes. 
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CHAPTER 3 

EXPERIMENTAL SETUP 

The experimental setup, as shown in Figure 4 in the form of a block 
diagram, consists of the following systems for a combustion jet ignition 
study 

1. Combustion Jet Igniter 

2. Main Chamber 

3. Flow Mixing System 

4. Optical System 

5. Electric Spark System 

6. Trigger Line 

3. 1 Combustion Jet Igniter 

The combustion jet igniter, as shown in Figure 5, was designed to gen- 
erate jets of combustion products containing free radicals and discharge 
them into the combustible medium as ignition sources. The igniter has a 
cylindrical chamber, so-called prechamber, of 1.59 cm in diameter and 
4.3? cm in length. The volume of the prechamber is 6.43 cm and can be 
varied by adding insert sleeves. The prechamber is fitted with gas inlet, 
purge outlet and jet orifice. The mixtures, coming from the flow mixing 
device, flow into the prechamber through the gas inlet located at the top 
of the prechamber. The mixtures can be purged from the prechamber 
through the purge outlet located at the side wail near the orifice plate. 
After continuous purging of the prechamber for about 3 minutes, the 
fresh mixtures are ignited, with both gas inlet and purge outlet valves 
closed. Ignition in the prechamber is achieved using an electric spark. 
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The electric spark occurs near the orifice plate between the valve rod 
head and the electrode protruding from the side wall of the inserted 
sleeve. The high electric voltage is applied through the high voltage lead 
to the valve rod. The valve rod head sits on the teflon seat housed in the 
orifice plate. The orifice plate housing the teflon seat can be inter- 
changed to provide variable-sized orifice diameters. 

The motion of the valve rod is controlled by a fast acting 19-SD-40 
solenoid. To close the valve, the valve rod is pushed down to seal the 
orifice by spring force. To open the valve, the valve rod is quickly pulled 
up by magnetic force induced by the solenoid. The motion of the valve 
rod is measured by a Bently Nevada displacement transducer of which a 
calibration curve is shown in Figure 6. 

Before the electric spark occurs, the gas inlet and the purge outlet 
valves are closed with a pressure tight seal, while the orifice is sealed by 
the valve rod head which is pushed down by a spring. Once the mixtures 
in the prechamber are ignited and combusted, the prechamber is pres- 
surized by combustion products. The pressure of the prechamber is 
measured by a Kistler Model # 601 pressure transducer and recorded on 
a Tektronix Model # 549 storage oscilloscope. The completely combusted 
products can be discharged into the main chamber through the orifice in 
the form of jets. This is done by activating the solenoid to open quickly 
the valve after the peak pressure is reached. 

By properly adjusting the timing for opening the valve before 
combustion is completed, i,e., before peak pressure is reached, the 
combustion jet igniter can generate jets of flames which are different 
from jets oi combustion products. In order to be able to study the 
difference between combustion jet ignition and flame jet ignition, a "spik- 
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ing circuit", as shown in Figure 7, was designed. This allowed the valve to 
be opened sooner enabling the incomplete combustion products to be 
discharged into the main chamber if desired. This spiking circuit powers 
the solenoid by initially raising the voltage across the solenoid coils to 
150 volts to overcome the coil's impedence. Once the spiking circuit is 
triggered, the current flows charge the capacitor and cut the high vol- 
tage, leaving 30 volts to keep the solenoid energized. This spiking circuit 
reduces the activation time, the time required for the solenoid to 
opening the valve, to 1 ms and cut the total opening time to 6 ms. 

3.2 Main Chamber 

The main chamber is a cylindrical, stainless steel explosion vessel 
fitted with optical glass windows 9 cm in diameter with a distance of 9 cm 
between their inner faces. The volume of the explosion vessel is 572 cm . 
and its viewing cross-sectional area is 64 cm The explosion vessel is 
fitted with four instrument plugs, one of which is at the top of the vessel 
and fitted with the combustion jet igniter as shown in Figure 8. Two addi- 
tional plugs v/ere used for purging and filling the vessel. The last plug was 
fitted with a Kistler Model #601 pressure transducer, the pressure signal 
from which was recorded on a Tektronix Model # 549 storage oscilloscope. 
The jets of combustion products, coming from the orifice of the 
prechamber, enter the main chamber fitted with optical glass windows 
suitable for Schlieren photography. The formation processes of combus- 
tion jets in the main chamber, filled with air at atmospheric pressure and 
temperature, were recorded on high-speed motion picture films. More- 
over, if the main chamber is filled with a combustible medium, then the 
ignition and inflammation processes of the medium caused by combus- 
tion jets can be investigated. 
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3.3 Flow Mixing System 

In order to measure ihe compositions of mixtures used in the 
prechamber, a continuous now mixing device was used for the sake of 
safety and easy operation. A schematic diagram of this system is shown 
in Figure 9. This system is an modified version of the continuous flow 
mixing apparatus designed by Kamel , which used hypodermic needles 
to measure flow rates. The advanced system uses jewel orifices to meas- 
ure How rate^ The flow is insured to be choked flow at the jewel orifice 
exit by maintaining a stagnation pressure upstream of the orii^e at least 
twice as large as the orifice downstream pressure. Under the critical flow 
conditions, the mass flow rate through the jewel orifice is ? linear func- 
tion of the upstream stagnation pressure and independent of the orifice 
downstream pressure. This linear dependence of critical mass flow rate 
on the orifice upstream stagnation pressure was calibrated for each jewel 
orifice using a "Precision" Wet Test Meter. The calibration curves for H ? , 
Og and C^Hq gases are shown respectively in Figures 10, 11 and 12. 

In an attempt to obtain better accuracy of mixture compositions, the 
continuous flow mixing system with 3 orifices was calibrated again using 
HP 5750B Gas Chromatograph. The calibration curves are shown in Fig- 
ures 13, 14 15 and 16. In obtaining these curves, the upstream stagnation 
pressures of any two orifices were maintained constant but that of the 
third orifice was allowed to vary. The lowest upstream stagnation pres- 
sure of these three orifices was verified as being at least twice the down- 
stream pressure of the orifices, thus insuring critical flow conditions for 
all orifices. 
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3.4 Optical System 

Schlieren photography was used in all combustion jet ignition experi- 
ments to record the formation and the inflammation processes of 
combustion jets. The main chamber, i.e. explosion vessel, was placed in 
the collimated beam of light created by a conventional 2 configuration 
Schlieren optical system, as shown in Figure 17. The mirrors used were 
152 mm in diameter and of 1.22 meter focal length. At the focus of the 
second n^rror was a circular aperture stop or a knife edge to discrim- 
mate the deflected light. Additional lenses and filters, not shown in Fig- 
ure 17, were used to insure proper imaging and exposure of the turbulent 
plume in the explosion vessel. The Schlieren movies were taken using a 
1000 watt Xenon continuous light source and a Hycam high speed camera 
filming at 5500 frames per second. The timing marks were made on the 
film by the camera every 1 ms. 

0.0 Electric Spark System 

The ignition circuit, as shown in Figure 18, used to energize the 
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combustion jet igniter, is similar to the one used by Morganstern et a L . 
The device uses a transformer to increase the voltage from 110 volts A. C. 
to 700 volts A. C.. Then, a rectifier converts the signal to 500 volts D. C. 
for the primary loop of a standard auto ignition coil. The magnetic field 
generated by the D. C. voltage is then eliminated by cutting off the 
applied vokage with an S.C.R. transistor switch. This generates up to 
40000 volts in the secondary loop of the ignition coil that can short- 
circuit the isolated spark gap and produce a spark in the prechamber. 
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3.6 Trigger Line 

As shown in Figure 4, a trigger signal enters the Hycam camera to 
activate the motion of the film. When the camera attains its pre-set 
speed for film motion, it sends out 2 signals to trigger respectively the 
oscilloscope and the ignition circuit. Then, the ignition circuit produces a 
spark in the prechamber and the oscilloscope gets ready for receiving 
the pressure transducer signals. In addition, the built-in delay circuit in 
the oscilloscope sends out a delayed signal to trigger the spiking circuit 
to activate the solenoid to open the valve The mixtures, with known 
compositions, are ignited by the electric spark and combusted to pres- 
surize the prechamber. The valve remains closed until the peak pressure 
is reached at which time it opens. Once the valve is open, the jets of 
combustion products are generated and enter the main chamber 
equipped with optical glass windows for taking Schlieren movies. The for- 
mation and the inflammation processes are recorded on the Schlieren 
movie film by the camera. In addition, the displacement-time trace for 
valve motion and the pressure-time traces for the prechamber and the 
main chamber are recorded on the oscilloscope. 
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CHAPTER 4 

EXPERIMENTAL RESULTS 

A series of tests were conducted to study the formation process of 
combustion jets. The experimental data obtained in this study consists of 
high speed Schlieren movies of combustion jet formation in the main 
chamber with synchronized pressure and valve motion records. Of par- 
ticular interest to us are the data of initial jet penetration speed which 
are listed in Table 1 These initial jet penetration speed were obtained 
from the linear regression analysis for the jet penetration depth data in 
the initial stage. The linear regression function for each test run is listed 
in Table 2. All of the results are grouped into 3 cases and are presented 
in the following sections. 

4. 1 Case I 

This case used a CgHg/Air mixture of equivalence ratio <p - 1.0 in the 
prechamber to study the formation processes of combustion jets. The 

g 

prechamber volume was 4.0 cm and its orifice diameter was 3 mm. 

Figure 19-a shows the pressure-time trace of combustion in the 
prechamber with all valves closed. The peak pressure for this test run is 
19 psig. This pressure-time trace is used as a reference for comparison 
in this case study. 

Figure 19-b presents the presure-time trace of combustion in the 
prechamber with valve opening starting at peak pressure. The spiking 
circuit shown in Figure 7 was not used in this test run. The 
displacement-time trace, which is an indication of valve motion, is aslo 
shown in Figure 19-b. The valve starts to open 4 ms after the electric 
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spark occurs and takes 7 ms to fully open. After the valve opens, the 
pressure decays at a rate faster than that of the reference pressure in 
Figure 19-a. The formation processes of the combustion jets for this test 
run are shown in Figure 20. The jets enter the main chamber, which is 
filled with air at atmospheric pressure and temperature, through the 3 
mm orifice in the form of a turbulent plume. The turbulent plume 
develops as it travels downstream. There are no flame fronts or 
inflammation processes observed in this turbulent plume. This is 
because the turbulent plume consists only of completely combusted pro- 
ducts. This is a typical case of combustion jets because the valve opens 
at or after the peak pressure. 

4.2 Case II 

In this case study, the spiking circuit was applied to reduce both the 
delay time required to start opening the valve and the valve opening 
duration. As shown in Figure 21 -a, the prechamber was pressurized up to 
24 psig with air, then the valve is suddenly opened using the solenoid 
applied with or without spiking circuit. The delay time required to start 
opening the valve is reduced from 4 ms to 2 ms and the valve opening 
duration was also cut from 7 ms to 4 ms. The prechamber pressure 
decays faster for the test run with spiking circuit than for that without 
the spiking circuit. This indicates that the air was discharged into the 
main chamber through the valve faster for the test run with spiking cir- 
cuit. These reductions in the time periods for the discharge of the 
combustion jets into the main chamber are very important in the 
combustion jet ignition study. 

This case again used a CgH 0 /Air mixture of tp = 1.0 to study the for- 
mation process of the combustion jet. The valve opening time is varied in 
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this case to investigate its effect on the turbulent plume formation 
processes and the jet penetration depth trace. As shown in Figure 21-b, 
the reference pressure-time trace for this case has a peak pressure of 19 
psig. All the results, as shown if Figures 22-27 for this case study are 
presented in the form of pressure and valve motion records with jet 
penetration depth trace. 

Tig are 22 shows the experimental results for a test run with valve 
opening 3 ms after peak pressure is reached. This time delay for starting 
valve opening after the peak pressure was found to be consistent witn the 
time period on the 16 mm movie film from the beginning time marker to 
the first frame of the jets. The pressure-time trace of this test rim repro- 
duced the reference pressure-time trace until the valve started opening. 
When the valve starts opening, the pressure is 11 psig and starts decayin? 
with the production of the combustion jet. The time period for generat- 
ing the combustion jet for this test run was 4 ms. The jet penetration 
depth data were obtained by measuring the length from the orifice to the 
turbulent plume tip which were recorded on the 16 mm Schlieren movie 
film. The regression line of the jet penetration depth data shows a pene- 
tration speed of 24.4 m/sec for the turbulent plume in the initial stage. 
But, when the turbulent plume travels closer to the bottom wall which is 
9 cm from the orifice, its penetration speed decreases mono tonic ally. 

Figures 23 and 24 show the experimental results for two test rims 
with valve opening 2 ms and 1 ms respectively after peak pressure. The 
pressures at the moment th valve starts opening for these two test runs 
are 15 psig and 17 psig respectively which are higher than last test run 
shown in Figure 22. The time period of generating combustion jets for 
these two test runs is 5 ms. The initial penetration speeds of the tur- 
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bulent plume for these two test runs are 29.5 m/sec and 30.6 m/sec 
respectively, and decrease monotonically when they travel closer to the 
wall. 

Figure 25 shows the experimental results for a test run with the valve 
opemng exactly when peak pressure is reached. This peak pressure of 19 
psig is the maximum pressure one can use to produce combustion jets 
for this case study. For this test run the time period for generating 
combustion jets is 6 ms and the intial penetration speed of the turbulent 
plume is 30.0 m/sec. The penetration speed again decreases when the 
turbulent plume travels closer to the wall. The jet penetration depth 
trace of this test run is compared with that of the test run shown in Fig- 
ure 22. This comparison, as shown in Figure 26, demonstrated that this 
test run with higher initial stagnation pressure upstream of the orifice 
has higher initial penetration speed of combustion jets. 

Figure 27 shows the experimental results for a test run with valve 
opening 2 ms before the reference peak pressure is reached. Thic is a 
test run with valve opening starting before combustion in the pechamber 
is completed. Therefore, combustion in the prechamber continues while 
the valve is opening. Finally, a peak pressure of 19 psig can be reached. 
The pressure at the moment the valve starts opening is 14 psig. The time 
period of generating combustion jets is 5 ms and the initial penetration 
speed of combustion jets is 27.7 m/sec for this test run. As the turbulent 
plume travels closer to the wall, its penetration speed decreases mono- 
tonically at a rate slower than all the previous test runs. It is interesting 
to note that the jet penetration depth trace for this test run is quite simi- 
lar to that for last test run shown in Figure 25. This indicates that contin- 
ued pressurization due to combustion in the prechamber after the valve 
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opens can maintain the initial jet penetration speed until the turbulent 
plume travels very close to the wall. Moreover, it is also found that flame 
fronts exist in the turbulent plume recorded on the 16 mm Schlieren 
movie film for this test run. Therefore, this test run should be con- 
sidered as a test run of flame jets. 

The results of this case study can be summarized as follows. The 
penetration speed of combustion jets is initially constant and decreases 
monotonically as the combustion jets travel closer to the wall. The con- 
stant penetration speed of combustion jets is proportional to the initial 
stagnation presure of the orifice for the same stoichiometric mixture. A 
flame jet can be generated with valve opening starting before the refer- 
ence peak pressure is reached. 

4.3 Case III 

This case used a CH^/Air mixture to study the formation process of 
the combustion. The equivalence ratio of the CH^/Air mixture is varied in 
this case to investigate its effect on the turbulent plume formation pro- 
cess. The spiking circuit was also applied in this case. 

Figures 28-a and 28-b show the reference pressure-time traces for 
the CH^/Air mixture of equivalence ratios <p = 1.0 and p = 1.2 respec 
tively. The peak pressure for <p - 1.0 is 18 psig and for p = 1.2 is 11 psig. 

The experimental data including the cinematographic Schlieren 
records of combustion jets for - 1.0 are shown in Figure 29 and 30. This 
test run of f - 1.0 starts valve opening at the peak pressue of 18 psig. 
The time period of generating combustion jets is 5 ms. The initial pene- 
tration speed of combustion jets is 22.7 m/sec and decreases monotoni- 
cally as they travel closer to the wall. The penetration depth trace of this 
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test run is compared with that of the test run shown in Figure 25. This 
comparison, as shown in Figure 31, demonstrates that combustion using 
a stoichiometric CgHg/Air mixture in the prechamber shows much higher 
initial penetration speed of combustion jets than that using a 
stoichiometric CH^/Air mixture although the initial stagnation pressure 
difference between these two test runs is 1 psi. This interesting com- 
parison needs further investigation. 

The experimental data including the cinematographic Schlieren 
records of combustion jets for <p - 1.2 are shown in Figures 32 and 33. 
This test run of ? - 1.2 starts the valve opening at the peak pressure of 
11 psig. The time period of generating combustion jets is 5 ms. The ini- 
tial penetration speed of combustion jets is 31.6 m/sec and decreases 
monotonically but not until they are much closer to the wall than in pre- 
vions cases. The penetration depth trace of this test run is compared 
with that of the test run shown in Figure 30. This comparison, as shown 
in Figure 34, demonstrates that combustion using a rich mixture in the 
prechamber tends to have higher initial penetration speed of combustion 
jets than that using a stoichiometric mixture although the former has 
lower initial stagnation pressure than the latter. This interesting obser- 
vation needs further study. 

The experimentaal data including the cinematographic Schlieren 
records of combustion jets for <p = 1.0 are shown in Figures 35 and 36. 
This test run starts opening the valve 1 ms after the electric spark. The 
combustion in the prechamber begins when the valve starts to open. 
Hence, the combustion cannot pressurize the prechamber up to a higher 
peak pressure. The peak pressure for this test run is 3 psig and the pres- 
sure at the moment the valve starts opening is 1 psig. The time period of 


generating combustion jets is 10 ms. The initial penetration speed is 19.5 
m/sec which is lower than any previous test runs. This initial penetration 
speed decreases monotonically when the turbulent plane travels closer to 
the wall. This is a test run of flame jets. Its jet penetration depth trace is 
compared with that of the test run shown in Figure 30. The comparison, 
as shown in Figure 37, indicates that although the initial stagnation pres- 
sure upstream of the orifice for flame jets is very low, flame jets have a 
comparable jet penetration speed as combustion jets. This interesting 
observation needs further study. 
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CHAPTER5 

THEORETICAL MODEL AND NUMERICAL TECHNIQUE 

In order to interpret and understand the physical phenomena 
observed in the turbulent jet experiments and furthermore to explain the 
formation processes of the large-scale eddy structure of turbulent jets, 
the theoretical studies of the formation and inflammation processes of 
turbulent jets were performed using Chorin’s Random Vortex Method 
implemented with the flame propagation algorithm. In this chapter, the 
theoretical models for the problems of turbulent jets with and without 
combustion are presented first, followed by descriptions of the numerical 
techniques used to solve the problems. 

5. 1 Theoretical Model 

The problem of turbulent jets in this study is simplified by assuming: 

1. The fluid is incompressible and two-dimensional; 

2. The pressure of the entire flow field is constant; 

3. Gravitational field is neglected; and 

4. The decoupled energy equation is neglected. 

With consideration of exothermic effects of combustion incorporated 
mto the turbulent flow field, two more assumptions are further made as 
follows: 

5. The flame with a specified constant normal burning relative velo- 
city is an interface where chemical reactions occur infinitely fast 
between reactants and products having constant but different 
densities; and 
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6. The exothermic effects of combustion on the turbulent flow field 
are manifested by the dynamic effect of combustion generated 
volumetric sources, flashing over the flow field, and accompanied 
with the transformation of reactants into products. 

After stating the above idealizations, the non-dimensional forms of 
mathematical models for the two cases without and with combustion are 
described in sections 5.1.1 and 5.1.2 respectively. 


5. 1. 1 Turbulent Jets without Combustion 


The non-dimensional governing equations, Navier-Stokes equations, 
for the 2-D, incompressible, free jet without combustion are: 

Continuity equation V u = 0 (5. 1) 

Momentum equation ~~= -VP + (5.2) 

L/t A| 


The initial and boundary 
specified as follows: 

Initial conditions 
Boundary conditions 


conditions for the problem are physically 

u s (o.O) for the entire flow field 
c = (0,0) at the solid boundaries (5.3) 
c = (0,1) at the inlet 


The problem with elliptical governing equations is well posed after all 
the peripheral boundary conditions and the initial conditions are 
specified. 


If vorticity u = Vxc is used, then the mathematical model is 
transformed into the following governing equation: 


Do 

Dt 




(5.4) 
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with D = 0 for the initial condition 

u = finite at all the boundaries satisfying the boundary 
conditions listed in Equation (5.3) (5.5) 

The problem of turbulent jets without combustion remaining to be 
solved is the governing equation (5.4) expressed in terms of u with the ini- 
tial conditions and boundary conditions listed in equation (5.5). The vor- 
ticity distributions, determining the turbulent jet flow field, at any time 
instant are unknowns in this mathematical model and were approximated 
by using Random Vortex Method The procedures of obtaining vorticity 
distributions at any time instant by using Random Vortex Method will be 
discussed in section 5.2. 

5. 1.2 Turbulent Jets with Combustion 

The mathematical model in sectibn 5.1.1 is further extended to 
include the combustion processes. Us*ng assumptions 4 and 5 mentioned 
in section 5.1, the complicated combustion processes are idealized as the 
transformation of reactants into products at the flame fronts with the 
combustion generated volumetric sources flashing over the whole flow 
field. Therefore, the exothermic effects of combustion are expressed in 
terms of the dynamic effects of combustion generated volumetric 
sources located at the flame fronts, la order to obtain the volumetric 
source distributions at any time instant, one has to know the locations of 
flame fronts which are described by the following flame propagation equa- 
tion: 

~ Sv ff/ (5-6) 

Now the mathematical model for turbulent jets with combustion can 
be set up by adding equation (5.6) into the mathematical model of section 
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5.1.1. The unknowns for the mathematical model with combustion are 
instantaneous flame front locations, i.e. volumetric source distributions, 
and instantaneous vorticity distributions which determine the turbulent 
jet flow field with combustion. The volumetric source distributions and 
vorticity distributions at any time instant in the flow field were approxi- 
mated by using Random Vortex Method implemented with the flame pro- 
pagation algorithm. The procedures of obtaining these distributions at 
any time instant will be discussed in section 5.2 

5.2 Numerical Techniques 

This section explains the numerical techniques and procedures used 
to approximate the solutions for the mathematical models listed in sec- 
tion 5.1. The section begins with the explanation of the Random Walk Con- 
cept and then performs the Schwarz-ChristofTel Transformation which is 
needed to consider the invisicid interaction of the vortices. Finally, the 
section explains how the Random Vortex Method and the flame propaga- 
tion algorithm were used to obtain the vorticity distributions and the 
volumetric source distributions at any time instant. 

5.2.1 Random Walk Concept 

The vorticity transport equation (5.4) indicates that if one follows the 
fluid element with a certain amount of vorticity, then he finds that the 
time rate of vorticity change is due to the vorticity diffusion. But in 
terms of the statistical point of view, Chorin proposed the Random Vortex 
Concept 30 ’ 0 to enable us to explain the equation (5.4) as the following 
idea. If one follows the fluid element with that specified amount of vorti- 
city, then he finds that the motion of the fluid element with a specified 
amount of vorticity is random. Here the fluid element of finite core and 
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with a small amount of vorticity is called a vortex blob and its diffusion 
processes are explained in terms of random motion based on a Gaussian 
distribution. The mathematical background for the idea of random 
motion, the so-called random walk, is discussed as follows: 

Consider initially imposing a fluid element with a specified amount of 
vorticity in an infinite motionless domain, then the vorticity transport 
equation is: 




(5.7) 


81 R, 

The solution for equation (5.7) with initial conditions u = <5 (r), a delta 
function at the origin, and boundary conditions u = 0 at infinity is 


, . , ,4rrl v _ z , R * rZ s 

u{r.t) = {—) exp ( — ) 


1 




r 1 / r 




1 _ / 1 r 

VSJf a 6Xp ” 2 a 2 


(5.0) 


This solution is a zero mean Gaussian Distribution with one standard 
deviation a ■VF . This distribution can be followed by the normal 
random variable, named the random walk of the vortex blob with vorticity 


u. The normal random variable, i.e. random walk, is further assumed to 
be time-independent to facilitate the study of transient formation 
processes of turbulent jets. Therefore, the random walk, which assumes 
the Gaussian Distribution of zero mean and one standard deviation 


a = -\/ -p- , at any time instant t = mbt can be obtained by m times the 

V R t 

summation of random walks which follow the Gaussian distribution of 
zero mean and one standard deviation a = Vtt 
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With the use of the random walk concept mentioned above, the dis- 
placement of a vortex blob was obtained at any time instant by adding 
the normal random displacement due to random walk to the expected 
displacement due to the inviscid interactions between vortex blobs. The 
procedures in obtaining the vorticity distributions at any time instant are 
discussed in section 5.2.3. 


5.2.2 Schwarz-Christoffel Transformation 

In order to facilitate the consideration of the inviscid interactions 
between vortex blobs or a vortex blob and a source, the Z-plane (physical 
plane) was transformed into the f-plane by performing the following 
Schwarz-Christoftel transformation: 

dZ _ 1 

where the Z-plane is non-dimensionalized by H and the £-plane by a. The 
Z-plane and f-plane were plotted in Figure 38 with points in the Z-plane 
being conformally mapped into those in the {--plane. As a result of this 
transformation, the uniform flow coming out of the inlet channel was 
mapped into a source located at the origin in the {--plane. 


In the {"-plane, the complex velocities for the following two- 

*S7 

dimensional simple flows are modified to be : 


_ i_ J_ 

d< ~ if < 




dW b 

d( ~ 2nMax(\f-( b \,r 0 ) (<•-<-*) 


Sc I <~<c 


1 


dW r 


dt; 2vMax{\$-t; c j,r„) 


for the source s located at the 
origin 

for any vortex blob 6 generated 
from walls 

for any volumetric source c gen- 


erated from combustion 
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where cj 6 is vortex strength and s c is volumetric source strength. The 

volumetric source strength s c , is determined^ 7 to be ' 

2 Pi at 

The complex velocities here are non-dimensionalized by the uniform velo- 
city C/ 0 the inlet of the channel. 

To determine the velocity of any vortex blob m induced by other vor- 
tex blobs, the source, combustion generated volumetric sources and 
their images due to the existence of the walls, the complex velocities for 
every simple flow are superimposed as follows: 


"SL. = U -iV = ^-+ V — 


* . - dW. 






where bg stands for the image of vortex blob 6 and eg stands for the 
image of volumetric source c. 


The induced velocity = U m +'iV TO ) of any vortex blob m is the 

velocity with which we follow to observe the random walk of this vortex 
blob m. This material velocity contributes to a displacement, so-called 
expected displacement in the statistical point of view, which is equal to 
dW 

times Af. The total displacement during the time period A* is 

approximated by adding a normal random displacement to this expected 
displacement. 


5.2.3 Random Vortex Method 

Tne Random Vortex Method has been successfully used to study the 
turbulent flows of different configurations. The detailed descriptions of 
this method can be found in references 53 and 57. Presented here are 
only the brief discussions of the use of this method to study turbolent 
jets. 
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As pointed out in section 5.2.1, .ie equation (5.4) statistically 
describes the motion of any wall generated vortex blob in the entire tur- 
bulent jet flow field. The motion of any vortex blob can be decomposed 
into two steps. First, the vortex blob undergoes an expected motion due 
to the inviscid interaction as discussed in section 5.2.2. Then, it has a 
random motion as described in section 5.2.1. While the expected displace- 
ment is calculated by the forward finite difference method as mentioned 
in section 5.2.2, the random displacement is determined by a randomly 
generated number from the computer. 

In order to update the vorticity distribution at any time step, the 
equation (5.4) is applied to every vortex blob in the turbulent jet flow field 
to describe its motion. After the vorticity distribution is updated, satis- 
faction of the no-slip boundary conditions with generation of vorticity at 
the walls is checked for the next time step. The generated vorticity is 
quantized and deposited at the walls This newly generated quantized 
vorticity becomes vortex blob once it move into the flow field. The vorti- 
city distribution with the newly generated quantized vorticity at the walls 
is updated again with the method mentioned above. These procedures of 
both checking the no-slip boundary conditions with generation of quan- 
tized vorticity and updating the vorticity distributions are repeated to 
obtain the vorticity distributions at any time step for turbulent jets. 

This Random Vortex Method uses the Prandtl boundary layer assump- 
tions to facilitate the calculations for the vortex blob interactions near 
the solid boundaries. The flow fields near the solid boundaries are called 
shear layer regions. The assumptions used in these shear layer regions 

are (1) —•« a nd (2) ^r« Hence the properties of a vortex 

53 

blob in these regions are simplified and the term vortex sheet is 
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applied instead of vortex blob. The newly generated quantized vorticity is 
called vortex sheet and is allowed to move in the shear layer regions. The 
vortex sheet diffuses or randomly walks only in the direction perpendicu- 
lar to the solid boundaries. In addition, the vortex sheet only influences 
other vortex sheets which are in the shadow region below it. These two 
properties of vortex sheets reduce a significant amount of computation 
time. Once vortex sheets move beyond the shear layer regions, they 
become vortex blobs. Conversely, vortex blobs will become vortex sheets 
if they move into the shear layer regions. A shear layer region thickness 
of 3 standard deviations was selected to give best accuracy in the calcula- 
tions. 

5.2.4 Flame Propagation Algorithm 

c Q 

The flame propagation algorithm was developed by Chorin 00 and suc- 
cessfully applied, with the use of the Random Vortex Method, to study the 
turbulent flow in a combustion tunnel . The algorithm was designed 
mainly to trace the motion of the flame fronts. Similar to equation (5.4), 
equation (5.6) describes the motion of the flame fronts in the turbulent 
jet flow field. The motion of a flame front can be decomposed into two 
steps. First, the flame front moves in the direction of its normal with a 
relative velocity taken here as a constant, equal to the appropriate lam- 
inar burning velocity of the mixture, S u . Secondly, the flame front 
moves with a velocity due to the inviscid ir*eractions as discussed in sec- 
tion 5.2.2. 

The motion of the flame front in the first step is due to consumption 
of the reactants and is called combustion. Combustion transforms reac- 
tants into products and generates volumetric sources located at the 
flame fronts. The combustion generated volumetric sources from the 
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previous time step are applied to influence the turbulent jet flow field at 
the current time step and then removed. Here, the exothermic effects of 
combustion on the turbulent jet flow field are manifested by the dynamic 
effects of combustion generated volumetric sources flashing over the 
entire flow field. The method of obtaining the combustion generated 
volumetric sources is incorporated in the procedures used in section 
5.2.3 to update the vorticity distribution at any time step for turbulent 
jets with combustion. 

The flame irum is an interface where chemical reactions occur 
infinitely fast between reactants and products. In order to locate the 
flame fronts, the flow field is divided into square cells with length Ac. The 
fraction of volume occupied by products in a given cell is defined as fol- 
lows: 

/ = '»rrvr = K‘ whereos / si 

Thus, / = 0 or 1 means that there are, respectively, either reactants 
or products in the cell, while 0</<l applies to cells containing the flame 
front interface. Its particular geometry is deduced, depending on the / 
numbers in neighboring cells. A detailed explanation of the interface 
geometry can be found in references 57 and 58. 
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-HAPTER 6 

THEORETICAL RESULTS 

With the use of the Random Vortex Method implemented with the 
flame propagation algorithm, four cases of computer simulations for tur- 
bulent jets were conducted to study the formation and the inflammation 
processes of turbulent jets. The typical input data used at the start of 
calculations are listed in Table 3. The calculated results of these four 
cases are discussed in the following sections. The purposes of studying 
these four cases are explained respectively in the following: 

1. Case I was designed to study the formation processes of the 
developing turbulent jets and the behavior of the large-scale 
eddy structure of turbulent jets at early time. This case simu- 
lates the formation processes of a turbulent plume created by 
the combustion jets, which are turbulent jets of combustion pro- 
ducts, before ignition. 

2. Case II was designed to study the behavior of the large-scale 
eddy structure of turbulent jets at later time. This case simu- 
lates the turbulent flow field of the developed combustion jets 
before ignition. 

3. Case III was designed mainly to study the inflammation processes 
in the turbulent jets. This case illustrates the mechanism of the 
inflammation processes in turbulent jets with a single point igni- 
tion source located at the center of the jet exit. 

4. Case [V was designed to study and simulate the inflammation 
processes initiated by a line source located at the jet exit. 
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6. 1 Case I 

This case was designed to study the formation processes of the tur- 
bulent plume observed in the experiments as discussed in Chapter 4. The 
calculated results are shown in Figures 40. Instantaneous vorticity distri- 
butions are shown on the left and instantaneous velocity vector fields are 
shown at six downstream locations of 0.5, 3.0, 5.5, 8.0, 10.5 and 13.0 from 
the jet exit. Nondimensionalized locations are based on the jet exit 
width. 

At the beginning of the computation (Time = 0.00), the entire tur- 
bulent jet flow is potential flow, which satisfies the initial condition of zero 
vorticity in the mathematical model. In the Z-plane, this potential flow is 
formed by discharging the uniform flow through the inlet channel into the 
semi-infinite plane as shown in the left of Figure 39. But in the 
corresponding <f-plane, the potential flow is formed by assuming a source 
located in the semi-infinite plane as shown in the right of Figure 39. 

Then as the calculations begin, the vortex sheets are generated by 
satisfying the non-slip boundary conditions at the solid walls. These vor- 
tex sheets are transported into the potential flow field as so-called vortex 
blobs. Vortex blobs of the same sense form a vortex cloud which has 
characteristics of self-rotating and expanding, causing the flow to recir- 
culate The vortex blobs produced from the horizontal walls are in the 
s^me sense as those produced from the adjacent vertical walls. The 
former vortex blobs were entrained into the main stream of the turbulent 
jets and combined with the latter vortex blobs to form two self-rotating 
and expanding vortex clouds. These two developing vortex clouds, which 
are of opposite sense, move parallel to each other and form the first vor- 
tex pair. 
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As the calculations continue, the vortex clouds continue growing 
separately by coalescing with the newly generated vortex blobs of the 
same sense and hence augmenting the entrainment. These coalescing 
processes are shown in Figures 40-a,b with this first vortex pair moving 
downward and its vortex clouds becoming bigger and bigger. Finally, the 
random nature of the RVM causes the first vortex pair to become anti- 
symmetric. While the first anti-symmetric vortex pair is growing down- 
stream of turbulent jets, the two incoming streams of vortex blobs are 
forced to move off the jet centerline and hence have time to form 
another anti-symmetric vortex pair before becoming entrained into the 
first vortex pair. This anti-symmetric effect of the vortex pair on forming 
another anti-symmetric vortex pair causes the vortex pairs to be stag- 
gered and forms a large-scale eddy structure of turbulent jets. 

It is interesting to note that the first vortex pair triggers the forma- 
tion of the large-scale eddy structure of turbulent jets. In this study of 
turbulent jets, the large-scale eddy structure of turbulence, as shown in 
Figures 40-c,d,e, is constructed by a series of vortex pairs, which are 
organized in the form of a staggered array of vortex clouds generating 
local recirculation flow patterns. 

This large-scale eddy structure of turbulent jets dominates the for- 
mation processes of a turbulent plume. All the fluids coming from the 
exit of the inlet channel are caused to recirculate by the large-scale eddy 
structure. The fluids from the left part of the orifice recirculate in a 
clockwise sense and the other part in a counterclockwise sense. More- 
over. the fluids from the region near the ihLt channel walls recirculate 
faster than those from the region far from the inlet channel walls. There- 
fore, the fluids exiting from the inlet channel are kept together under 
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recirculation and form a turbulent plume. The penetration speed of this 
turbulent plume is determined by the strength of the first vortex pair 
heading the large-scale eddy structure of turbulent jets. 

In this case calculations were performed only up to 21 time 
steps. This case study has demonstrated not only the formation processes 
of the first vortex pair and turbulent plume, but also the mechanism of 
developing a large-scale eddy structure. The continuation of calculations 
becomes very time-consuming and costly because of the large number of 
vortex blobs (about 4600) produced. In order to continue the calcula- 
tions and to study the behavior of the large-scale eddy structure of tur- 
bulent jets at later times, one has to reduce the generation rate of vortex 
blobs. The effects of the non-slip horizontal walls on the turbulent jet 
flow fields are already known from this case study to enhance the 
development of the first vortex pair with a larger degree of randomness, 
and hence to augment the large-scale eddy structure of turbulent jets. 
The assumption of frictionless horizontal walls will be used in the studies 
of the following cases. 

6.2 Case II 

In designing this case, the horizontal walls, which were non-slip in the 
previous case, were assumed to be frictionless here. Hence, the genera- 
tion of the vortex blobs from the horizontal walls was eliminated to facili- 
tate the caculations. In an attempt to obtain a steady flow, the calcula- 
tions for this case were conducted up to 66 time steps to simulate the 
turbulent flow field created by the developed combustion jets. Tb calcu- 
lated results are shown in Figures 41 on a scale of 6/10 to those of Fig- 
ures 40. Instantaneous vorticity distributions are shown on the left. 
Instantaneous velocity vector fields are shown on the right at six 
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downstream locations of 0.0, 4.0, 0.0, 12.0, 16.0, 20.0 from the jet exit. 
Nondimensionalized locations are based on the jet exit width. 

Figures 41-a,b,c,d show the developing processes of the first vortex 
pair as explained in case I. The first vortex pair in this case is formed by 
two clouds, which are produced only from the vertical walls, of less ran- 
domness and smaller size. Therefore, the potential core region of the tur- 
bulent jets is longer in this case, and the large-scale eddy structure of 
the turbulent jets is formed further downstream 

As the calculations continue, the two big vortex clouds, which are 
always associated with the existence of the non-slip solid boundaries, 
form an anti-sy:mmetric first vortex pair and start introducing the second 
smaller vortex pair. Gradually, as shown in Figures 41-e,f,g, the second 
vortex pair develops l. to an anti-symmetric form due to the anti- 
symmetric effect of the first vortex pair. The second vortex pair also 
introduces the third vortex pair, while itself being entrained into the first 
vortex pair. The second vortex pair finally coalesces with the first vortex 
pair to form a bigger, more anti-symmetric first vortex pair. The third 
vortex pair hence becomes the second vortex pair and is influenced by 
the first vortex pair to become more anti-symmetric. As calculations 
continue further, the third, the fourth, and subsequent vortex pairs are 
formed, as shown in Figures 41-h,i,j,k,l,m,n, using the same mechanism 
as the second vortex pair. This series of vortex pairs is becoming bigger, 
more anti-symmetric due to the anti-symmetric effect resulting from the 
preceding vortex pair and by coalescing with the successive vortex pairs. 
The vortex pair formed near the jet exit finally will be entrained into the 
first vortex pair by coalescing itself with the preceding vortex pair. This 
process enables the development of the large-scale eddy structure, which 
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is organized in the form of a staggered array of vortex clouds, of greater 
amplitude further downstream. This large-scale eddy structure of tur- 
bulent jets dominates not only the entrainment and mixing processes but 
also the entire turbulent flow field. 

6.3 Case III 

This case was designed to study the mechanism of the inflammation 
processes in turublent jets. A single point ignition source was applied at 
the center of the jet exit at the beginning of the calculations. The results 
of calculations are shown in Figures 42 with the instantaneous vorticity 
distributions on the left and the instantaneous flame front locations on 
the right. The scale of these figures is same as that of Figures 41. 

As the calculations begin, the large-scale eddy structure develops in 
a mechanism similar to that described in Case II. In addition, the forma- 
tion of the large-scale ec’dy structure in this case is under the influence 
of the exothermic effects of combustion. The large-scale eddy structure 
for this case is augmented in such a way that the first vortex pair is 
bigger in size, and the potential core region is shorter in length. Hence, 
the inflammation processes can enhance the formation processes of the 
large-scale eddy structure of turbulent jets. 

As discussed in Cases I and II, the large-scale eddy structure dom- 
inates the entire turbulent jet flow field. In this case, the large-scale 
eddy structure also demonstrates its control over the flame propagation 
processes. That is, the flame fronts, with a prescribed normal burning 
speed (S u = 0.02), are carried by the turbulent flow field. Of particular 
interests is that the flame fronts are carried by the recirculatiing flow 
fields created by the large-scale eddy structure of turbulent jets. 
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Besides, the recirculation flow fields also keep the burnt mixture 
together to form a turbulent plume with combustion. The penetration 
speed of this turbulent plume is determined, instead of by the normal 
burning rate, primarily by the strength of the first vortex pair heading 
the large-scale eddy structure of turbulent jets enclosed within the flame 
fronts. Therefore, the inflammation processes and the flame shape for 
this case of relatively small normal burning speed are determined by the 
large-scale eddy structure, as well as where and when the ignition source 
occurs in the turbulent flow field. 

This case performed calculations up to 21 time steps and illustrated 
the effects of inflammation on turbulence, and vice versa. For the sake of 
brevity, the case only used a single point ignition source located at the 
center of the jet exit to study the mechanisms of the inflammation 
processes in turbulent jets and the formation processes of a turbulent 
plume with combustion. The next case considf . o 5 point ignition sources 
to study the inflammation processes initiated by combustion jets. The 5 
point ignition sources form a line ignition source at the jet exit. 

6.4 Case IV 

This case was designed with an attempt in mind to simulate the 
inflammation processes initiated by combustion jets. The case assumes 
that the mixtures coming out of the exit are complete combustion pro- 
ducts. Moreover, the case assumes that no induction period is needed for 
the initiation and propagation of the self-sustained flames. In addition, 
this case assumes to have a line ignition source at the exit of turbulent 
jets. The calculations begin with the flames starting to propagate from 
the jet exit. The results of calculations are shown in Figures 43 with the 
instantaneous flame fronts superimposed on the instantaneous vorticity 
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distributions. The scale of these figures is the same as that of Figures 41. 

As the calculations continue, the large-scale eddy structure develops 
by the same mechanism as described in Case III. The turbulent plume 
with combustion also develops in the same way as described in Case III. 
The large-scale eddy structure again can dominate the entire turbulent 
jet flow field and the inflammation processes in it. 

As the calculations continue further, the large-scale eddy structure 
develops itself into a staggered array of vortex clouds which create local 
recirculation flows, causing the turbulent plume to have more wrinkled 
flame fronts. The recirculation flows also keep the burnt gases together, 
as they did in Case III, but producing a flame shape which is different 
from that in Case III. This illustrates that the flame shape and 
inflammation processes have a strong dependence on the number and the 
locations of the ignition sources deposited into turbulent jets. The pene- 
tration speed of a turbulent plume with combustion is also determined, 
instead of by the normal burning rate, primarily by the strength of the 
first vortex pair heading the large-scale eddy structure of turbulent jets 
enclosed within the flame fronts. 

The case performed calculations up to 36 time steps and again illus- 
trated the effects of the large-scale eddy structure on the flame shape 
and inflammation processes of turbulent jets. 
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CHAPTER 7 
CONCLUSIONS 

The conclusions for this combustion ,jt igniU i study are summar- 
ized below: 

1 . A combustion jet ignition system was developed to generate tur- 
bulent jets of combustion products containing free radicals and 
tc discharge them as ignition sources into the combustible 
medium. 

2. Experimentally, the penetration speed of combustion jets has 
been found to be constant initially and decreases monotonically 
as turbulent jets of combustion products travel closer to the 
wall. 

3. The initial penetration speed of combustion jets is proportional 
to the initial stagnation pressure upstream of the orifice for the 
same stoichiometric mixture. 

4. Combustion using a CgHg/Air mixture in the prechamber shows 
higher initial penetration speed of combustion jets than that 
using a CH 4 /Air mixture although the initial stagnation pressure 
difference between these two test runs is 1 psi. This is a primi- 
tive study and needs further investigation. 

5. For combustion using a CH^/Air mixture in the prechamber, the 
initial penetration speed cf combustion jets is higher for rich 
mixture than for stoichiometric mixture although the former has 
lower initial stagnation pressure than the latter. This is a primi- 
tive study and needs further investigation. 
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6. The Random Vortex Method is a very useful tool for studying the 
large-scale eddy structure of turbulent jets. 

7. The first vortex pair triggers th? formation of the large-scale 
eddy structure of turbulent jets. 

8. The large-scale eddy structure of turbulent jets is constructed 
by a series of vortex pairs, which are organized in the form of a 
staggered array of vortex clouds generating local recirculation 
flow patterns. 

9. The large-scale eddy structure dominates the formation 
processes of turbulent jets or experimentally observed turbulent 
plumes. 

10. The flame fronts are carried by the recirculating flow fields 
created by the large-scale eddy structure. In addition, the recir- 
culation flow fields keep the burnt mixture together, forming a 
turbulent plume with combustion. 

11. The large-scale eddy structure dominates the inflammation 
processes in turbulent jets. 

12. The flame shape and inflammation processes have a strong 
dependence on the number and the locations of the ignition 
sources deposited into turbulent jets. 
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Table l 

Jet Penetration Speed for Combustion Jets Experiments 


Case 

Figure 

Mixture 

<P 

ti (ms) 

Pw(psig) 

K(m/sec) 

II 

22 

CgHg/Air 

1.0 

+3 

11 

24.4 

II 

23 

CgHg/Air 

1.0 

+2 

15 

29.5 

II 

24 

CgHg/Air 

1.0 

+ 1 

17 

30.6 

II 

25 

CgHg/Air 

1.0 

0 

19 

30.8 

II 

27 

CgHg/Air 

1.0 

-2 

14 

27.7 

III 

30 

CH 4 /Air 

1.0 

0 

18 

22.7 

III 

33 

CH 4 /Air 

1.2 

0 

11 

31.6 

III 

36 

CH 4 /Air 

1.0 

♦ 

1 

19.5 


cp - Equivalence ratio. 

Po, - Initial stagnation pressure upstream of the orifice. 

V i - Initial penetration speed. 

- Valve opening time after (+) or before (-) peak pressure. 
♦ - Valve opening 1 ms after electric spark. 
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Table 2 

Linear Regression Functions 
(Y = AT + B) 


Case 

Figure 

A(m/sec) 

B(m) 

r 

II 

22 

24.38 

2.53xl0" 3 

0.996 

II 

23 

29.47 

-9.20xl0' 5 

0.996 

II 

24 

30.62 

2.05x1 0‘ 4 

0.996 

II 

25 

30.00 

3.26xl0* 4 

0.994 

II 

27 

27.67 

2. 02x1 O' 3 

0.999 

III 

30 

22.69 

3.02X10' 3 . 

0.997 

III 

33 

31.56 

1.14xl0‘ 3 

0.999 

III 

36 

19.50 

- 9.49xl0* 4 

0.998 


r - Correlation coefficient 
Y - The estimated jet penetration depth 

A - Slope of the least squares line of the jet penetration depth data 

B - Y-intercept of the least squares line of the jet penetration depth 
data 

T Time (second) 



Table 3 


Input Data for Turbulent Jets Calculations 


Reynolds number 

10000 

Width of jet exit 

1.0 

Time increment 

0.1 

Space increment 

0.2 

Length of each vertical wall 

0.2 

Length of each horizontal wall 

4.5 

Side width of a cell 

0.2 

Normal burning speed 

0.02 

Density ratio (Reactants /Products) 

6.0 

Strength of a vortex blob 

0.05 

Minimum strength of a combustion 


generated source 

0.001 

Standard deviation 

4.47 x 

Shear layer thickness 

1.34 x 
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Fig. 1 
Fig. 2 


Fig. 3 

Fig. 4 
Fig. 5 
Fig 6 
Fig. 7 
Fig. 0 
Fig. 9 
Fig. 10 
Fig. 11 
Fig. 12 
Fig. 13 
Fig. 14 
Fig . 1 j 


FIGURE CAPTIONS 

Idealized Otto cycle 

Comparison of pressure-time traces 

Curve 1 - Obtained in closed bomb experiments using conven- 
tional spark ignition 

Curve 2 - Obtained in closed bomb experiments using combus- 
tion jet ignition 

Curve 3 - Obtained from theoretical calculation for constant 
volume combustion using CEC-76 computer program 

Jet fluid field in the plane of symmetry of the jet. Reynolds 
number = 2500, d = 0.3 in, 0< x/d <35 (from Dimotakis et al ) 

Block diagram of experimentaLsetup 

Schematic diagram of combustion jet igniter 

Calibration curve for a displacement transducer 

Schematic diagram of the spiking circuit 

Main chamber fitted with combustion jet igniter 

Schematic diagram of the continuous flow mixing device 

Critical flow rate calibration curves for Ng-orifice 

Critical flow rate calibration curves for Og-orifice 

Critical flow rate calibration curves for CgHg-orifice 

Calibration curve of CgHg/Og for continuous flow mixing device 

Calibration curves of Ng/Og for continuous flow mixing device 

Calibration curves of CgHg/Air for continuous flow mixing devjce 
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Fig. 16 
Fig 17 
Fig. 18 
Fig. 19- 

Fig. 19- 

Fig 20 

Fig 21- 
Fig 21- 

F ig 22 


Calibration curves of CH^/Air for continuous flow mixing device 

Schematic diagram of Schlieren optical system 

Schematic diagram of ignition circuit for combustion jet igniter 

Pressure-time trace for combustion of CgHg/Air mixture (? = 
10) in the prechamber with all valves closed (Reference 
pressure-time trace for case I) 

Pressure-time and displacement-time traces for combustion of 
Cg Hg/Air mixture ($f = 1.0) in the prechamber with valve open- 
ing at peak pressure (without spiking circuit) 

Cinematographic Schlieren records of combustion jets in the 
main chamber filled with air at atmospheric pressure and tem- 
perature (without spiking circuit) 

Mixture : CgHg/Air (? = 1.0) 

Orifice diameter : 3 mm 

The effect of the spiking circuit on the valve motion and the 
pressue-time trace 

Pressure-time trace for combustion of CgHg/Air mixture (^ = 
1.0) in the prechamber with all valves closed. (Reference 
pressure-time trace for case II) 

Experimental data for a test run with valve opening 3 ms after 
peak pressure 


Mixture : C^Hg/Air ( tp = 1.0) 
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Fig. 23 Experimental data for a test run with valve opening 2 ms after 

peak pressure 
Mixture : CgHg/Air (? = 1.0) 

Fig 24 Experimental data for a test run with valve opening 1 ms after 
peak pressure 

Mixture : C^Hg/Air (? = 1.0) 

Fig 25 Experimental data for a test run with valve opening at peak 
pressure 

Mixture : CgHg/Air {<p - 1.0) 

Fig. 26 Comparison of jet penetration depth traces between two test 
runs with valve opening at peak pressure and 3 ms after peak 
pressure 

Fig. 27 Experimental data for a test run with valve opening 2 ms before 
peak pressure 

Mixture : CgHg/rtir = 1.0) 

Fig. 28 Pressure-time traces for combustion of Cri^/Air mixture in the 
prechamber with all valves closed (Reference pressure-time 
traces for case III) (a) <p = 1.0 (b) <p = 1.2 

Fig. 29 Cinematographic Schlieren records of combustion jets in the 
main chamber filled with air at atmospheric pressure and tem- 
perature (with spiking circuit) 

Mixture : CH^/Air (? = 1.0) 

Orifice diameter : 3 mm 


Valve opening : at peak pressure 
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Fig 30 Experimental data for a test run with valve opening at peak 
pressure 

Mixture : CH^/Air (<p = 1.0) 

Fig. 31 Comparison of jet penetration depth traces between two test 
runs with CH^/Air mixtures and CgHg/Air mixtures 

Fig. 32 Cinematographic Schlieren records of combustion jets in the 
main chamber filled with air at atmospheric pressure and tem- 
perature (with spiking circuit) 

Mixture : CH^/Air (<p = 1.2) 

Orifice diameter : 3 mm 
Valve opening : at peak pressure 

Fig. 33 Experimental data for a test run with valve opening at peak 
pressure 

Mixture : CH 4 /Air = 1.2) 

Fig. 34 Comparison of jet penetration depth traces between two test 
runs with CH^/Air mixtures of equivalence ratios ? = 1.0 and = 
1.2 

Fig. 35 Cinematographic Schlieren records of combustion jets in the 

main chamber filled with air at atmospheric pressure and tem- 
perature (with spiking circuit) 

Mixture : CH^/Air (<p = 1.0) 

Orifice diameter : 3 mm 

Valve opening : 1 ms after electric spark 
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Fig 36 Experimental data for a test run with valve opening 1 ms after 
electric spark 
Mixture : CH 4 /Air {<? = 1.0) 

Fig 37 Comparison of jet penetration depth traces between two test 
runs with valve opening at peak pressure and 1 ms after electric 
spark 

Fig. 38 ochwarz-ChristofTei transformation of Z-plane into f-plane 
Left : Physical piane or Z-plane 
Right : Transformed plane or f-plane 

Fig 39 Streamline patterns of initial flow of turbulent jets in Z-plane 
(left) and f-plane (right) 

Fig 40 Computer plots ot turbulent jets without combustion for R t - 
10000 

Left : Instantaneous vorticity distribution 

Right : Instantaneous velocity field 

(a) Time = 1 vb) Time = 6 (c) Time =11 (d) Time = 16 

(e) Time ~ 21 

Fig. 41 Computer plots of turbulent jets without combustion for R, = 
10000, assuming smooth horizontal walls 

Left : Instantaneous vorticity distribution 

Right : Instantaneous velocity field 

(a) Time = 1 (b) Time = 6 (c) Time =11 (d) Time = 16 

(e) Time = 21 (f) Time = 26 (g) Time = 31 (h) Time = 36 

(i) Time =41 (j) Time = 46 (k) Time - 51 (1) Time = 56 


✓ 


& 


r 



63 


(m) Time = 61 (n) Time = 66 

Fig. 42 Computer plots of turbulent jets with combustion for R t = 

10000, Su = 0.02, smooth horizontal walls and a single point igni- 
tion source located at the center of the jet exit 

Left : Instantaneous vorticity distribution 

Right : Instantaneous flame front locations 

(a) Time = 1 (b) Time = 6 (c) Time =11 (d) Time = 16 

(e) Time = 21 

Fig. 43 Computer plots of turbulent jets with combustion for R, = 
10000, S u = 0.02 smooth horizontal walls and a line ignition 
source located at the jet exit 

(a) Time = 1 (b) Time = 6 (c) Time = 11 (d) Time = 16 
(e) Time = 21 (f) Time = 26 (g) Time = 31 (h) Time = 36 



64 


oRvaw^ 1 - 


OF POOR Q' 



Process 1-2 : 
Process 2-3 : 
Process 3-4 : 
Process 4-5 : 

Process 5-6 : 

Process 6-1 : 


Isentropic compression stroke 
Constant- volume combustion 
Isentropic expansion stroke 
Constant-volume and adiabatic exhaust 
blowdown 

Constant-pressure and adiabatic exhaust 
stroke 

Constant-pressure and adiabatic intake 
stroke 


Figure 1 
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Figure 5 
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RV - Regulating Valve 
MV - Manual Valve 
SV - Solenoid Valve 


Figure 9 
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JO - Jewel Orifice 
P - Pressure Gauge 
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ORIGINAL PAGE IS 
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0.15 mm 
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for C^Hg-orifice 
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Figure 12 
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Figure 13 



Stagnation pressure upstream of 0 5 -orifice 
is fixed (=57 psig) * 



Figure 1 4 





Figure 15 
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ORIGINAL PAGE B 
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As a 



Figure 17 

















ORIGINAL PAGE W 
OF POOR QUALITY 



Displacement 


Pressure 


a 


Time Scale : 2 ms/div. 
Pressure Scale : 5 psi/div. 



Dj splacement 


Pressure 


b 


Time Scale : 2 ms/div. 

Pressure Scale : S psi/div. 
Displacement Scale : 1.13 mm/div. 
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a 

Time Scale : 2 ms/div. 

Pressure Scale : 10 psi/div. 
Displacement Scale : 1.13 mm/div. 



b 

Time Scale : 2 ms/v.Iv. 
Pressure Scale : 5 psi/div. 


Pressure 

Displacement 


Displacement 


Pressure 


Figure 21 
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Displacemeu 


Pressure 


Time Scale : 2 ms/div. 

Pressure Scale : 5 psi/div. 
Displacement Scale . 1.13 mm/div. 





Jet Penetration Depth (cm) 


ORIGINAL PAGEJ3 
OF POOR QUALITY 



Time Scale : 2 ms/div. 

Pressure Scale : 5 psi/dlv. 
Displacement Scale : 1.13 mm/div. 






Jet Penetration Depth (cm) 




Jet Penetration Depth (cm) 






o- Valve opening at peak pressure 

•- Valve opening 3 ms after peak pressure 
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Figure 26 
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Figure 27 
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Displacement 


Pressure 


a 

Time Scale : 2 ms/div. 
Pressure Scale : 5 psi/div. 
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Time Scale: 5 ms/div. 
Pressure Scale : 5 psi/div. 


Figure 28 
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Figure 29 
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Time (ms) 


o- C^Hg/Air mixture 
• - CH , /Air mixture 


Figure 31 


• o 










Jet Penetration Depth (cm) 


ORIGINAL PAGE la 
OF POOR QUALITY 
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Time Scale : 5 ms/div. 

Pressure Scale : 5 psi/div. 
Dispalcement Scale : 1.13 mm/div. 
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Figure 33 
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Figure 36 
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Figure 37 
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